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Abstract 

"“^Research  completed  Includes  the  following: 

The  dynamics  of  a  system  of  nonlinear  ordinary  differential 
equation  depends  on  the  constant  coefficients  (parameters)  of  the 
system.  Identifying  these  coefficients  from  the  solution  curves 
defines  an  Inverse  problem.  A  method  to  determine  the  values  of 
the  parameters  from  a  finite  number  of  solution  curves  was 
developed  and  Implemented.  The  method  consist  of  two  major 
algorithmic  procedures' (1)  A  derivative  free  nonlinear 
optimization;  (2)  An  error  analysis  of  the  parameters  found. 


The  Optimization  Algorithm  [ 1 ] 

The  nonlinear  optimization  algorthm  utilizes  random  vector 
as  directions  of  search  to  find  an  optlmun  point.  The  components 
of  the  random  vector  are  independently  generated  from  a  Gaussian 
distribution.  Two  Interpolating  schemes  are  used;  (a) 
lagranglan  polynomial  approximations,  (b)  and  spline 
approximations.  The  method  is  iterative  and  has  fast  convergence 
when  used  on  problem  having  many  variables  (more  than  ten)  .  The 
algoirthm  is  capable  of  finding  an  optimun  point  even  when  the 
function  F  to  be  optimized  is  not  available  in  closed  form,  but 
rather  only  values  of  F  at  discrete  points  can  be  obtained. 
Moreover,  derivatives  of  F  are  not  available.  This  kind  of 
functions  are  typical  when  dealing  with  inverse  problems  where  a 
set  of  parameters  has  to  be  determine  from  a  finite  number  of 
discrete  points  on  the  solution  curves.  [  2  ] 
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The  method  does  not  require  "close'*  estimates  of  the  optlmun 
point,  and  It  is  easy  to  implement  and  use.  The  algorithm  was 
tested  on  several  dynamical  systems  and  nonlinear  functions  in 
many  variables,  such  as  a  model  of  gluconeogener is  having  31 
parameters,  [3].  Also,  the  "Rosenbrootc"  function  of  50  variables 
1 4 J  and  the  "Powell  singular"  function  [5]  which  has  the 
characteristic  that  precisely  at  its  minimum  value  the  Jacobian 
becomes  singular  (this  is  the  reason  that  Gradient  methods  fall 
to  converge  to  the  minimum).  The  results  on  the  test  problems 
showed  the  versatility  of  the  method,  and  its  superior 
performance  compare  to  often  algorithms  [1]. 


The  Error  Analysis  [ 6 ) 

An  important  aspect  in  the  parameter  estimation  technique  is 
the  validation  of  the  parameters  found  by  the  optimization 
technique.  Small  perturbations  in  the  observations 

r  “  l,...,m  can  result  in  a  percentage  error  that  can  vary 
greatly  (by  orders  of  magnitude)  between  parameters.  Thus  1 
performed  an  error  analysis  of  the  parameter  values  to  determine 
the  validity  of  the  results.  The  methods  can  be  outline  as 
follows:  Denote  by  K  e  R^  the  best  estimate  found  by 

optimizing  the  function  F,  let  {X]  -  [f(X,K,t)  ,[X(0)]  -  [C] 

S  S  S  8 

describe  the  system  for  the  initial  condition,  and  let 


i 
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(D(t)l  -  [A(t)J  l D( t ) ]  +  [ B ( t ) )  ,  [0(0) ]  -  0 

“  8  8  B  S 

be  the  corresponding  variational  system,  where 

We  integrate  the  variational  system  at  t^,  ...  ta  and  form 

the  matrix  H  such  that 

"  '  .-1  Li  • 

where  Wr  is  a  weighting  function  for  each  data  point. 

Then  ,  the  expected  variance  for  the  i1*1  parameter,  is 

i 

given  by 


r!!il 

LiSkjJs  * 


[D( t) ] 


*lr  (Hii)_l 


To  minimize  the  probability  of  making  a  mistake  during  the 
derivation  of  [D(t)]g,  I  implemented  an  algorithm  that  uses 
automated  symbol  manipulation  to  formally  obtain  the 
[bf  ^/hkj,  )]  g,  [  (bf  ^/bx ]  ,  and  the  necessary  sum  and  product  of 

such  matrices  [6]. 
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